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1. Introduction 

Global expectation values of the chiral condensate and quark number density give information 
about the structure of the thermal medium at finite temperature and non-zero chemical potential /i. 
Calculating these observables locally, i. e. depending on their spatial position, allows for a spatially 
resolved analysis of dynamical quark clouds and the chiral condensate around static quark, diquark 
or Meson sources. We investigate correlation functions of the chiral condensate/quark number 
density with the Polyakov loop or products of Polyakov loops, which allow to probe the screening 
and string breaking mechanism of systems composed of static color charges at finite temperature. 
Lattice calculations of these operators have been performed in the past with staggered and Wilson 
quarks [|l], ||, ^ for rather large quark masses and on small lattices. The investigations presented 
here form the basis for future studies of the spatial correlators of the quark number density, which 
has been proposed as a means to decide whether the QGP has a liquid phase [^]. 

2. Operators 

We have calculated the chiral condensate and the quark number density depending on their 
spatial position x by using 

ZZ^(x) = ^Tr, IM- and <(x) = ^Tr, IM^, (2.1) 

where x = {xq,x), f = l,...,nf labels the flavor, /i/ is the corresponding chemical potential and 
Trc denotes the trace in color space, with Trc 1 = 1. Derivatives with respect to /X/ are evaluated 
at iJ.f = 0. We can now calculate the correlator of XX^{^) «{(x) with the Polyakov loop L(x) 
or products of Polyakov loops, for which we use 

L(x)=Z(g2)^^Trc n ^4(xo,x), Ceg(r)=L(0)r(r), CggW = L(0)L(r), (2.2) 

^0=0 

where L, and therefore also Cqq and Cqq, is properly renormalised by the multiplicative renor- 
malisation constant Z{g^)^\ which is determined at T = from a matching of the static quark 
potential to the string potential. The Cqq and Cqq describe static meson and diquark systems, 

respectively. The correlation functions ^«{(x)L(y)^ and (nq{x)Cs{r)'^ for s = QQ,QQ are then 
properly renormalised as well. The renormalisation for XX^ involves both multiplicative and ad- 
ditive renormalisation constants. In order to eliminate the additive contributions, we calculate the 
connected correlator 

{XI^{x)L{y))^ = (zZ^(x)L(y)> - {xx^) (L) , (2.3) 

which approaches the corresponding continuum expression (yy-^ (x)^(y))^. in the limit a —^0. 

In this work, we analyze gauge field configurations generated in 2-1-1 flavor QCD with the 
p4fat3 staggered fermion action, where the strange quark mass is fixed to its physical value 
and the light quark mass is choosen as m/ = O-lm^. Lattice sizes are 16^ x4,6 For the 

calculation of expectation values we employed the random noise vector method with 100 random 
noise vectors. 
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Figure 1: {xX^^c (upper row) and (n^'^L) (lower row) for the strange quark species. 

3. Results for {xxL)c (ngL) 

In fig. |l] we show the resuhs for {xX^L,)c (upper row) and (m^Z.) (lower row) for the strange 
quark species obtained from Nt = 4 lattices. We observe that all correlators vanish rapidly at 
distances rT>\, i. e. when XX^ ^^'^ approach their thermal expectation value undisturbed by 
the presence of the static quark. For small separations {xX^^c assumes smaller values, signalling 
a gradual restoration of chiral symmetry in the vincinity of the static quark. The correlator {n'^Lj 
becomes smaller than zero for smaller separations as well, as the static quark source is screened by 
an excess of dynamical anti-quarks arranged in a cloud around the static quark. The effect becomes 
smaller above T^, as thermal gluons are available in the deconfined phase to provide screening. All 
these findings are in qualitative agreement with those of earlier studies [|l|, ^, ^, Q]. 

For {n^Lj we observe a peculiar dip at rT ^ 0.5, which we investigate in more detail in fig. 0. 
Varying the lattice spacing, we see in fig. ^ (left), that the position of the dip varies with a as well, 
pointing at a lattice artifact that vanishes in the continuum limit. In fig. |^ (right) we computed 
(?i*L) at T /Tc = 0.905, A^T = 4 for p4, Naik and naive staggered discretisations of the fermion 
matrix M. The artifact is present in all three discretisation schemes, thus its appearance is not 
dependent on the particular action used. Rather the mixing of different quantum channels in the 
staggered formulation is probably responsible for this effect. 

In order to quantify further the behavior of these correlators, we have performed simple expo- 
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Figure 2: for Nz — A and 6 (left) and for different staggered discretisations of the fermion matrix M 

at T/T, = 0.905, = 4 (right). 



nential fits. The resulting screening masses m''^ of the correlator involving the chiral condensate 
are, of course, not affected by the multiplicative renormalisation constants. In fig. || we show 
{XX^^c (upper row) and (n^qL) (lower row) for the strange quark species on a log-scale together 
with the results of the ;t^^-fits of exponential functions. We obtain similar results for the light quark 
species. The fit ranges employed were 0.3 <rT<\ for (riqL^, where we left out the values due 

to the artifact at rT ^ 0.5, and 0.1 < < 1.2 for (xx^L)^. 

The results for the screening masses obtained from the fits with exponential functions are 
shown in fig. ^left). As expected^ the screening masses obtained from (riqL^ are larger than those 

from {xX^^)c ^^'^ generally find is larger than m'. We observe a drop of all screening 
masses close to Tc, as all correlation functions become large in the cross over region. For large 
temperatures we have m''" /T ^ 1 and we therefore expect m' = to hold, which is not yet the 
case at T/r^ = 2. Moreover, the deviation of the data from the lattice with smaller aspect ratio 
points to volume effects which are still present in our calculations of m' *. The Debye masses 
extracted from the Nt = 6 lattices are roughly of the same size as the screening masses [^. In order 
to avoid contamination from excited states, we plan to study p = projections of the correlators in 
the future. 

The total induced quark number of flavor / around a static anti-quark can be computed by 

where we show the results of our calculation in fig. fright). As was already observed in 2-flavor 
QCD [lIo|], (Nq j_ vanishes rapidly for T > Tc, where screening is predominantly provided by ther- 



Q 

mal gluons. For T —fO the induced quark number approaches one, where this limit is approached 
faster for light quarks than the strange flavor, since the heavier sea quarks are more difficult to 
induce. 



'For heavy sea quarks the correlation functions {xX^ ^"9^^ '-^'^ related to hopping parameter expansions 
of real and imaginary parts of Polyakov loops, respectively!^]. In this limit one finds at high temperatures m„ lI^xxL = 
1.5 [il. 
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Figure 3: {xX^^c (upper row) and (nji) (lower row) for the strange quark species on a log-scale. Lines 
are the results of the ;)f^-fits with simple exponential functions. 



4. Results for (n^Cgg) and (n^Cgg) 

We now turn to the computation of three point correlation functions (n^Cgg) and (n^yCgg). 
In fig. Q we show (Mj^(xi,X2)CQQ(r)) for T /Tc = 0.905, calculated on lattices with temporal extent 
Nt = 4 for separations rT = 1.25 (left) and rT = 0.25 (right) of the Polyakov loops, which are 
placed at xi = ±rT /2,X2 = 0. We see that for large separations r the correlator assumes negative 
values centered around the individual positions of the static quarks. This signals an excess of 
dynamical anti-quarks and the screening mechanism is mesonic, i. e. QQ —>■ Qq + Qq, where Q 
stands for a static quark and q for a dynamical quark. For small separations we observe positive 
values of the correlator. Thus the screening mechanism has changed to baryonic, i. e. QQ QQq. 
This behavior has been observed in the induced quark number of the diquark system in 2-flavor 



QCD as well [|T§]. 

In fig. ^ we show the values of («'^Cgg) (left) and ("^Cgg) (right) along the line connecting 
the static quarks for several temperatures and two separations, where xi = x,X2 = 0. The dotted 
lines mark the position of the Polyakov loops, where the anti-quark source in the 22-system is 
located at larger values of x. We see mesonic screening in the diquark system for large separations 
r at all T calculated here, whereas at the smallest distance only the lowest temperature turns to 
baryonic screening. In the 22-system the screening mechanism is mesonic for all T and rT, 
i- e. QQ^Qq + Qq. 
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Figure 4: Left: Screening masses obtained from an exponential fit. Closed symbols: Nx — 4, open symbols: 
Nx = 6. Debye masses are taken from [^. Right: Induced quark number around a static anti-quark. 
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Figures: (n^Cgg) for T/r, 
Polyakov loops. 



: 0.905, A^T = 4 at separations rT = 1.25 (left) and rT = 0.25 (right) of the 



5. Conclusion and Outlook 

We have studied correlators of the local chiral condensate and the local quark number density 
with Polyakov loops in 2+1 flavor QCD with almost realistic quark masses at finite temperature on 
16^ X 4, 6 lattices. We found that the correlator of the Polyakov loop with the quark number density 
and the connected correlator with the chiral condensate are well described by an exponential fit- 
Ansatz. We calculated the screening masses from these correlators and the induced quark number 
in the static quark system. Moreover, we investigated the quark number density in static diquark 
and mesonic systems, where we found the screening mechanism in the diquark system to change 
from mesonic to baryonic at small separations for small temperatures. 

In the future we want to study the p = projection of the correlators with the Polyakov loop to 
avoid contamination from excited states in the screening masses. Furthermore we like to examine 
lattices with a larger physical volume in order to get better control over finite volume effects. This 
work forms the basis for future computations of spatial correlators of the quark number density 
(^nq{0)nq{r)^, which might shed light on whether the QGP has a liquid phase 
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Figure 6: Cut along the connecting axis of the static quarks, {n-l^Cgg) (left) and (ngCgg) (right) for several 
temperatures. The dotted Unes give the position of the Polyakov loops. 
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